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Abstract. We derive some rigorous results concerning the backflowr operator 
introduced by Bracken and Melloy. We show that it is linear bounded, self adjoint, 
and not compact. Thus the question is underlined whether the backflow constant 
is an eigenvalue of the backflow operator. From the position representation of the 
backflow operator we obtain a more efficient method to determine the backflow 
constant. Finally, detailed position probability flow properties of a numerical 
approximation to the (perhaps improper) wave function of maximal backflow are 
displayed. 



PACS numbers: 03.65 

1. Introduction and summary 

Let a 1-dimcnsional free solution of the Schrodingcr equation contain positive mo- 
menta only, and let Px(t) be this wave function's probability (at time t) to detect the 
particle at any position > x. Then Px{t) starts out from at time t = — oo and tends 
towards 1 for i — > oo. Because of Px{t) = j{t,x), the (position probability) current 
j(t,x) is naively expected to be nonnegative for every {t,x). Yet there exist positive 
momentum wave functions such that the current at, e.g., a; = is negative at certain 
intermediate times. In this case the half space probability as a function of time, i.e., 
Pq : R ^ [0, 1] is not monotonically increasing. 



This so called quantum backflow effect seems to have been mentioned first by 
AUcock in his work on the time of arrival in quantum physics ^ , while Bracken and 
Melloy [2] have given the first detailed account of the phenomenon in 1994. AUcock 
presented the backflow effect in order to disprove the hypothesis that the current at 
a; = yields the probability density of arrival times for a free positive momentum 
wave packet at a; = 0. Recently it has been shown that the backflow effect indicates 
discrepancies among two other proposals of arrival time densities. More specifically 
it has been shown in reference [2] that none of the arrival time densities, which obey 
Kijowski's axioms j^, coincides with the one of Bohmian mechanics [S]- Furthermore 
their average arrival times differ if and only if the wave function in question leads 
to backflow, in which latter case the average Bohmian arrival time precedes that of 
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Kijowski's distributions. 



Bracken and McUoy [21 posed the question whether the backflow of probability is 
restricted by a stronger bound than the obvious one given by 1. Though the existence 
of such a stronger bound was not to be expected, they attempted to numerically com- 
pute the smallest upper bound A for the decrease of P. By converting this backflow 
constant A into the supremum of the spectrum of an integral operator K in momen- 
tum space, surprisingly enough, Bracken and Melloy approximately found its value 
to be 0.04. Meanwhile the precision of the value of A has been improved by Eveson, 
Fewster, and Verch to 0.038452. 

In the present work we describe a new approximation method to determine A, 
which provides independent confirmation of the results of . Such confirmation is in 
need since a rigorous proof for the conjecture A < 1 is still missing. The basic idea is 
to use a decomposition of the integral operator K into a sum of Fourier transformed 
multiplication operators. In this way the method of fast Fourier transform becomes 
applicable and A can be approximated with less computational effort. We obtain an 
improved value for A of 0.0384517. As a byproduct of our numerical computations 
wc approximate the (perhaps improper) wave function of maximal backflow and we 
exhibit some of its more detailed position probability flow properties. 

The primary goal of this work, however, is to provide some exact results 
concerning the integral operator K of Bracken and Melloy. From a unitary equivalence 
it will become obvious that K is linear bounded and self adjoint. Then we prove that 
K is not compact by showing that —1 belongs to the spectrum of K yet it is not an 
eigenvalue. We have not been able to conclusively answer the question whether A is an 
eigenvalue of K in the strict mathematical sense. However we shall provide numerical 
plausibility that this is indeed the case. A more extensive discussion of some of our 
results concerning the backflow phenomenon is given in reference 0. 



2. The backflow constant 



The free Schrodinger evolution Ut : (M) (K) from time to time t e M is 

given in the momentum representation by 

{Ut^) (fc) = cf>t (fc) := exp {-ikh) <t>{k). 

Here t denotes the rescaled time variable htphys/{2m). Let tpt denote the inverse 
L^-Fourier transform of 0t, i.e., 

1 

ipt{x) (.F*0t) {x) = —= / exp (ifca;) (f>t{k)dk. 
V 27r J-oo 

Let a particle have the momentum space wave function (/) at time 0. If Ijc/ill = 1, the 
probability that a position measurement at time t yields a position a; > reads 



Here H : (R) — > (M) denotes the orthogonal projection with 

(n/)(x) = ( ^[f) f'^^^^ _ 

^ ■> ' 1 for x < 
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If a unit vector <j> E L'^ (E) has its support contained in IR>o, i.e., if 110 ~ 0. the 
probabihty P{(t)t), according to DoUard's lemma [Sj, obeys P{4>t) — > for t — > — oo 
and P{4't) ~* 1 for t — !■ oo. However, the mapping t ^ P{4't) need not monotonically 
increase from to 1. Rather it may decrease during several intermediate time 
intervals. [5] Thus there exist momentum space wave functions (f> G 7i+ := 11 (L^ (M)) 
such that P{4>s) > P{'Pt) holds for some s < t. For such holds 

A {(f>) := sup {P{(t>s) ~ P{(i>t) I s, i e R with s <t}> 0. 

Unit vectors (p E H+ without backflow yield A (</>) = 0. We define the backflow 
constant by 

A := sup{A((/)) I e with ||(/)|| = 1}. 

Introducing the orthogonal projection Ht TUT* Ut we obtain for any unit 

vector (M) 



p{(t,s) - P{(f>t) = ^<^, (n, - lit 

Because of 

- n, = uu. ( ffi^ - n^^l Ui±s 

-2- \ 2 2 J 2 

it follows that 

A = sup{^0,c/; (n^T-nr) c/r^) \ (f> en+,\\(t)\\ =i,reR,re 

Since the unitary Ut stabilizes ?i+ we infer 
A = sup 

T>0 

where Bt denotes the backflow operator 

Bt -^Il-T -^T, (1) 

and a (A) denotes the spectrum of a linear operator A. This follows from theorem 2, 
section 8, chapter XI of !) . Observe the bounds —id < Bt < id. 

Let the one parameter family of unitary dilation operators : (R) (M) 
with /i G R>o be given by (V^^) (fc) = y/jL(() {jik) . The operators commute both 
with n and with TTiT* and a brief computation shows that 

v^Utv; = u^2t. 

From this it follows that 

V^JIBtIW; = IiB^2TlV. 

Since the spectrum of an operator is invariant under a unitary transformation we have 
the following result, on which our numerical computation will be based. 

Proposition 1 For any fixed real T > holds A = super (HBtH) ■ 

In view of this result we choose T = f in what follows. The corresponding 
operators Ut=i and Bt=i will be abbreviated by U and B. 
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3. Equivalence with the treatment of Bracken and Melloy 

Now we will prove that our definition of A indeed is equivalent to the one of Bracken 
and Melloy These authors heuristically introduce A via time integrals of currents 
at point X = over arbitrary finite intervals. From this they motivate their final 
definition of A as the supremum of the spectrum of the integral operator 

K : (R>o) ^ (M>o) 

with (Kf) (fc) = -i r (,) d,. 



Jo k- q 

Let rj : H+ {^>o) denote the unitary operator with {ricf)) (fc) = (f){k) for all 

fc > and for all cj) in H+. 

Proposition 2 For all (f> G 7i+ there holds Krjcf) = j]IlBIl(p, i.e., the restriction of 
HBH to TL+ and K are unitary equivalent. 

Proof. Since liBIi is bounded it is sufficient to show ■t^IiBI{(t> = A'jyt/) for all from a 
dense subspace T) C Wc shall choose T) = 5+(M), the space of all C°° functions 
from R to C with fast decrease and with their support contained in M>o- 



1 f°° 

H : (R) ^ (R) , {Hf){k) = -V 

J-oo 



-dq. 



As a prerequisite we first demonstrate a relation between the orthogonal 
projection J-HJ-* and the Hilbcrt transformation 

k ^ q 

Here V indicates that the improper integral is meant as the principal value. For 
/ G 5(K.) wc obtain by means of Lebesgue's dominated convergence theorem and by 
means of Sochozki's formula ^01 

{TllF*f) {k) =^j^ (^1^ e'^V {q) dg) da: 

^-:^hmr^M_d. 
27ri £\o q - k + ie 



2m 



By continuity wc infer 



THJ'* = i {-iH + id) . (2) 

From equation ^ it is easy to show that the Hilbcrt transformation is unitary and 
that a{H) = {i, -i}. 

From the equations and ^ follows 

B = ^ {UHU* - U*HU) . (3) 
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From this we obtain for fc > and (f) G V 

inBIl4>) (fc) = ^ iHU*<j>) (fc) - ^ (HU^) (fc) 

1 /-^ sin (P-q^) , , _ , 

=--y^ ^--ii0(,)d,.(i.0)(fc). 

Clearly for fc < holds (UBIlcf)) (fc) = 0. By continuity we have HBYlcf) = 7]^^Kr](j) for 

all (/) e 7i+. Thus the restriction of HBH to Ti.+ is unitary equivalent to K. 

□ 

Therefore the defining relation of 0, A = supa {K) , indeed holds. 



4. Noncompactness 

Proposition 3 The backflow operator IlBIl is not compact. 

Proof. For HBU holds -id < UBU < id. Therefore a (UBU) C [-1,1]. For every 
unit vector </> S T-i+, according to Bollard's lemma [S] holds 

lim {(f>,IlBTllcf)) = -1. 

T — >oo 

Since the spectrum of HBt^ does not vary with T it follows that — 1 G a(nBir). 
If UBU were compact, then —1 were an eigenvalue of nSII. Let cf) £ L'^ (R) with 
= 1 denote an eigenvector of IlSn with eigenvalue —1, i.e., HBUcj) = —<f) holds. 
Since IIBIIc/) S Ti+ it holds that (f) G 7i+ . Then it follows from the triangle inequality, 
from the unitarity of the Hilbcrt transformation H , and from equation (PJ that 

1 = II0II = linsn^ii = lins^ii 

< \\B4>\\ = ^ \\{UHU* - U*HU) 011 

< i {\\UHU*q^\\ + \\U*HUcf>\\) = i {\\HU*^ + WHUcPW) = 1 
Thus the triangle inequality becomes an equality and we have 

\\{UHU* - U*HU)(j)\\ = \\UHU*(j)\\ + \\U*HU4>\\ 
from which it follows that there exists some a G C such that 
UHU*(t) = aU*HU(t). 

Since UHU* and U*HU are unitary it follows that \a\ — 1. From the above sequence 
of inequalities it also follows that 

||ni?0|| = IIB0II . 

This is equivalent to liBcf) = B(f>. Thus the eigenvector condition HBYlcj) = —(p implies 
B(j) = —(j), from which by means of equation (O it follows that 

Thus C/0 =: (f> G 7i+ is an eigenvector of H. Since ct {H) = {i, — i} it follows that 
a — 1 G {2, —2}. Because of |q;| = 1 this implies a = — 1 and if $ = i$. Thus it follows 
that 

TWT*<^ = i {-\H + id) $ = $. 
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Thus we have njF*cf> = for some nonzero $ G 7Y+. Now the following lemma 

implies the contradiction $ = 0. Thus —1 is not an eigenvalue of the backflow operator. 
Since every nonzero spectral value of a compact operator is an eigenvalue, the backflow 
operator necessarily is noncompact. □ 

Lemma 1 Let <i> e (R) with n<I> = $ and UT*<^> = T*<i>. Then $ = holds. 

Proof. Any function from (R) is locally intcgrablc. Therefore the inverse Fourier 
transform of $ e Ti.+ is the distributional boundary value of the holomorphic function 
$ on the complex upper half plane deflncd by 

~ 1 P°° 

$ : {z e C I 3z > 0}, a + ib^—= e'"'=e-^'=$(fc)dfc. 

V Stt Jo 

If the boundary value obeys ITF*$ = T*^, then the distribution JF*$ is zero on 
]R<o- From the generalized uniqueness theorem, see theorem B.IO on p. 100 of |12| . 
it follows that 4> = 0. Thus also the boundary value JF*<I> of $ vanishes. Since !F* is 
unitary we also have $ = 0. □ 



5. Numerical Computation of the backflow constant A 

In 121 E] the integral operator K is approximated by a finite square matrix, whose 
largest eigenvalue is taken as an approximation of A. If, however, we apply the power- 
method to the expression for A, which is given in proposition ^ we immediately 
approximate the largest eigenvalue without having to compute any matrix. One only 
needs to apply multiplication operators and fast Fourier transformations to an arbi- 
trary initial vector. The power method works as follows, jll) 

Let the matrix A e C^^^ be symmetric. Let a be the eigenvalue of A with the 
largest absolute value. Let Vq G be a nonzero vector with nonzero component 
within the eigenspace of A corresponding to a. Then the sequence ivn)n£Tig 
recursively defined by 

1 



\Vn\\- 



-Av„ 



Then holds 



1™ 4+1 • -ir~ii' 

n^oo \\v„\\ 



Since ct (IIBII) C [— 1,A] we apply the power method to the nonnegative, dis- 
cretized operator HBH + id. Its largest eigenvalue then approximates A + 1 while Vn 
tends towards the corresponding eigenvector. 

The analysis was started with A^o = 10'' grid-points covering the interval [0,q'o], 
where go is set to 50. Now the power- method was applied with 1000 iterations to a 
constant starting vector. Then we repeated the computation for to = A'o'^ grid- 
points and a larger momentum interval [0,q] with q = qoVh for h ~ 1,2, . . .40. In 
this way the covered interval grows while the absolute step size q/N gets smaller. The 
results Xh for different factors of accuracy h then were used to extrapolate to /i ^ oo 
leading to an approximation Aoo for the backflow-constant. The results of this com- 
putation can be seen in flgure^ 
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In order to check for the possibihty that the constant starting vector has van- 
ishing component within the eigenspace of the dominating eigenvalue various other 
starting vectors have been chosen as well. After only few iterations they all led to 
the same results. Since it seems extremely unlikely that all chosen starting vectors 
have vanishing components within the eigenspace of the dominating eigenvalue, our 
algorithm is likely to approximate the largest spectral value of HBH + id. 
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Figure 1: A plotted against h and fit Aoo + b/^/h. 
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Figure 2: A plotted against 1/y/h and polynomial fit of third order. 

An even better result for Aoo is achieved by fitting the graph of figure El with a 
polynomial of third order. The extrapolated value for the backflow-constant can then 
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be read off from the intersection of the y-axis with the graph. (This corresponds to 
the point 1/Vh = on the x-axis and h ^ oo respectively.) This yields 

Aoo = 0.0384517 

which agrees with the value given in |H] by 0.038452. 

By means of the power-method we also get an approximation of the possibly 
improper eigenvector associated with the backflow constant. It will be discussed briefly 
in the next section. 

6. Approximate backflow maximizing vector 

Since the operator K of Bracken and Melloy is real, the (improper?) backflow maxi- 
mizing eigenvector may be chosen to be real valued in the momentum representation. 
From this it follows that the position representation at time has even real part and 
and odd imaginary part. More generally, the time dependent wave function is invari- 
ant under the combined parity and time reversal operation. 

We take as an approximate backflow maximizing vector the vector w„ obtained 
from the power method, where we choose N = 10^, q = 50 and we make n = 1000 
iterations. The starting vector vq is - as before - simply the constant function. This 
leads - as one can read off from figure El ^ to quite a bad approximation of A by about 
0.0297, but a further increase of the accuracy leaves the appearance of the approxi- 
mate eigenvector i;„ as displayed in figure 13 totally unchanged. 



0.1 
0.05 





■-10 -8 -6 -4 -2 2 4 6 8 10 



0.02 
0.01 



-0.01 
-0.02 

-15 -10 -5 5 10 15 

Figure 3: v„ in momentum and configuration space (below, real part 
, imaginary part ). 

The position probability density of w„ subject to the free time evolution is dis- 
played in figures 0] |S1 These figures by themselves do not provide unquestionable 
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0.2 




Figure 4: Position probability density of for —20 < x < +20 at 

times t e {-3, -2,..., +3}. 




Figure 5: Position probability density of {un)t for —20 < a; < +20 and 
-3 < t < 3. 



In order to strikingly illustrate the backflow we compute for the approximate 
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backflow maximizing vector u„ the current j{t,x) at position x = as a function of 
time. The result is shown in figure El where the backflow-domain ] — 1, 1[ is plainly 
indentifiable. This interval seems to be the only one in which ij„ leads to a backflow 
through a; = 0. The area below it, as required, approximately sums up to the backflow- 
constant. The corresponding half space probability as a function of time is also shown 
in figure El Further evidence for the backfiow phenomenon of u„ is provided by figure 
d This figure shows some integral curves of the space time vector field (l,j/p), the 
flow lines of the Bohmian velocity field, within the backfiow-domain. All the integral 
curves which pass the line a; = at a time t with — 1 < t < 1 pass it in the negative 
direction. 
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Figure 6: j{t, 0) ( ) of Vn and the corresponding half space 

probability ( ) as functions of time. 

The question which remains open is this: Is there really a backfiow-eigenvalue 
- in the strict mathematical sense - to which Ago is an approximation? From the 
approximate eigenvector ii„ evidence can be found that there is indeed one. To this 
end we compute the contribution of the interval [0, g] to the norm-square of w„ and 
compare it to /g' |/(fc)pdfc with /(fc) = M ■ sin(fc^)/A: where A/" is a normahzation 
constant. Note that / € The results are shown in figure |S1 The two graphs 

are very similar and the norm of the w„ seems to converge even faster than that of /. 
Thus it seems plausible that A is indeed an eigenvalue of the backflow operator. 
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Figure 7: Flow lines of the velocity field j/p of The fiow lines 

between pass through x = in the negative direction. Two 

consecutive flow lines are separated by a probability of approximately 
2.4- 10-3. 
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Figure 8: Norm-squares of Vn ( ), and sin(fc^)/fe ( ). 
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